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We reply to the Comment of Brodsky and Gardner on our paper “New limits on intrinsic charm 
in the nucleon from global analysis of parton distributions” [Phys. Rev. Lett. 114, 082002 (2015)]. 

We address a number of incorrect claims made about our fitting methodology, and elaborate how 
global QCD analysis of all available high-energy data provides no evidence for a large intrinsic charm 
component of the nucleon. 


In a recent Comment [1], Brodsky and Gardner (BG) 
criticize our global PDF analysis [2] of all available high- 
energy scattering data, including those from fixed-target 
experiments at high x and low Q 2 , which placed strong 
constraints on the magnitude of intrinsic charm (IC) in 
the nucleon. For a range of models of IC, the analysis [2] 
strongly disfavored large magnitudes of IC, with the mo¬ 
mentum fraction carried by charm quarks (*} IC at most 
0.5% at the 4er confidence level (CL). 

BG claim that because our global analysis [2] uses 
0(30) parameters, as is typical in all such fits, one must 
adopt a much larger tolerance criterion than Ay 2 = 1. 
In fact, it is well known that for Gaussian distributions 
parameter errors in y 2 fits are determined by Ay 2 = 1, 
irrespective of the number of parameters in the fit [3, 4], 
The parameter m in Table 38.2 of Ref. [3], for example, 
is the dimensionality of the error regions for joint distri¬ 
butions (m = 1 for linear errors, m = 2 for error ellipses, 
etc.), and has nothing to do with the total number of 
parameters in the fit. Actually, Fig. 38.2 of Ref. [3] in¬ 
volves the number of degrees of freedom of a fit (number 
of points — number of parameters) and not the number 
of parameters in the fit. For the determination of indi¬ 
vidual parameter errors, the correct dimension is m = 1, 
which gives Ay 2 = 1 at the 68.3% CL. (For examples of 
error ellipses with m = 2, see Fig. 12 of Ref. [5].) 

The parameter errors and y 2 profiles related to one- 
dimensional probablility distributions are correctly eval¬ 
uated using Ay 2 = 1. Errors on other quantities are then 
computed using standard error propagation techniques, 
such as the Hessian method; they can also be used to pro¬ 
duce error regions of different dimensionalities with the 
appropriate Ay 2 criteria [3, 4]. Apparently, BG have con¬ 
fused the dimensionality of error regions with the number 
of independent parameters in a fit. Their claims about 
Ay 2 are simply wrong. 

Tolerance criteria Ay 2 > 1 are used by some PDF 
groups [6-8] on purely phenomenological grounds, to ac¬ 
count for tensions among different data sets, while others 
[5, 9, 10] use the standard Ay 2 = 1. The y 2 profiles in 
[2] were presented as a function of (*} IC , so that (*} IC 


values for different tolerance choices can be easily com¬ 
pared. BG also suggest that our single parameter errors 
were obtained by fixing the other parameters at the y 2 
minimum. This is not true: we minimize the y 2 with 
respect to all other parameters in the fit, as is standard 
procedure in global fits. Had we not properly refitted the 
complete model, the rise in the y 2 away from the mini¬ 
mum would be even steeper than for the profiles shown 
in Fig. 1 of Ref. [2], 

Inclusive DIS cross sections, such as those measured 
at SLAC, receive contributions from all quark flavors, so 
they cannot by themselves provide significant constraints 
on charm. The power of a global fit, however, lies in 
the correlation between different observables, with dif¬ 
ferent weightings of quark flavors, within the framework 
of perturbative QCD. While the bulk of the data from 
SLAC [11] at large x lie below the charm threshold, cross 
sections below threshold constrain light quark distribu¬ 
tions, which indirectly impacts on the determination of 
IC at the same kinematics. Our analysis also accounts 
for the suppression of charm production below and near 
the hadronic charm threshold [1, 2]. Implementing the 
suppression involves some model dependence in relating 
the partonic and hadronic charm thresholds [2, 6], and 
while this affects the quantitative limits (with partonic 
threshold factors alone ( x) lc would be < 0.1% at the 5cr 
CL), the effects do not alter the overall conclusions about 
the magnitude of IC supported by the data. 

To avoid dealing with complications from thresholds 
and other hadronic effects at low W 2 and Q 2 , many 
global PDF analyses impose more severe cuts on W 2 and 
Q 2 than those in Ref. [2], While this simplifies the theo¬ 
retical treatment, it also removes a significant amount of 
data at large x that could potentially impact the question 
of IC. 

Recently, some PDF analyses [5, 7, 10] have relaxed 
the conventionally more restrictive W 2 and Q 2 cuts in 
order to better constrain large-* PDFs. Such analyses 
benefit from increased statistics at large x, but require 
careful treatment of subleading 1/Q 2 and nuclear correc¬ 
tions. Our analysis [2] employs the standard treatment of 


2 


target mass corrections (TMCs) [12], phenomenological 
higher twists determined consistently within the same fit 
[5], and the latest technology in nuclear corrections [7]. 
Apparently confusing Refs. [5] and [7], BG assert that 
we model higher twists as isospin independent, and that 
our TMCs are problematic at x —> 1. In fact, our higher 
twist corrections do depend on isospin, as evident from 
Table III of Ref. [5] , and are determined empirically with¬ 
out assuming any functional form. 

Furthermore, the well-known threshold problem of 
TMCs at x = 1 is relevant only at very low W 2 , well 
below the cuts made in all global PDF analyses [12]. It 
is also not true that we neglect intrinsic strangeness and 
bottom: the s and s PDFs are parametrized model inde¬ 
pendently at the input scale, and given our results for IC, 
intrinsic bottom is negligible for the current phenomenol¬ 
ogy [13]. 

Our global fit carefully propagates all statistical and 
systematic errors, both uncorrelated and correlated, in¬ 
cluding normalization, for all data sets used. (For details 
of the fitting code see Ref. [ L4].) For the SLAC data, our 
analysis uses the original hydrogen and deuterium cross 
sections [11] rather than the derived structure functions 
(obtained by combining measurements at different ener¬ 
gies) , which allows for a more exact treatment of point to 
point correlated errors. Aside from the SLAC data, other 
measurements, such as the NMC proton and deuteron 
cross sections [15] and the inclusive proton cross sections 
from HERA [16], also disfavor nonzero values of IC. 

In addition to the fit of the standard high energy data 
sets used by most PDF groups, in Ref. [2] we also con¬ 
sidered a fit including data from the EMC measurement 
of the charm structure function F!f [17] — sometimes 
cited [18] as providing evidence for large IC in the nu¬ 
cleon. In practice, the EMC data have strong tension 
with other measurements, and give a very large overall 
X 2 /7Vciat of > 4, and a Q 2 dependence incompatible with 
perturbative QCD. Several of the EMC data points at the 
highest x values (x > 0.2), where there are no other di¬ 
rect constraints from charm production experiments, lie 
systematically above all global fits, including ones with 
IC contributions [2]. At the same time, at low x val¬ 
ues (x < 0.02) where charm distributions are strongly 
constrained by HERA [19], however, the EMC data are 
significantly below the fitted results. 

We thus disagree with the assertion of BG that (x) IC ~ 
0(1%) is “consistent with the analysis of the EMC mea¬ 
surements” [1]. No reasonable amount of nuclear correc¬ 
tions (which are, in fact, considered in Ref. [2]) or Ay 2 
tolerance can reconcile the EMC F% data with the rest 
of the global data set within a QCD framework, without 
invoking a very peculiar shape for IC that is strongly at 
variance with all IC models [2, 20, 24]. Consequently, no 


modern QCD analysis [5-10, 21-23] includes the EMC 
data in their fits. The MSTW analysis [6] compared F% 
computed from their PDFs with the EMC charm mea¬ 
surements, and concluded that “If the EMC data are to 
be believed, there is no room for a very sizable intrinsic 
charm contribution.” We agree with this conclusion. 

This material is based upon work supported by the 
National Science Foundation (T.J.H. and J.T.L.) under 
Grant PHY-1205019. The work of T.J.H. was also sup¬ 
ported in part by DOE Grant No. DE-FG02-87ER40365 
and DE-FG02-97ER-41014. P.J.-D. and W.M. were sup¬ 
ported by the DOE Contract No. DE-AC05-06OR23177, 
under which Jefferson Science Associates, LLC operates 
Jefferson Lab. 


[1] S. J. Brodsky and S. Gardner, arXiv:1504.00969. 

[2] P. Jimenez-Delgado, T. J. Hobbs, J. T. Londergan and 
W. Melnitchouk, Phys. Rev. Lett. 114, 082002 (2015). 

[3] K. A. Olive et al., Chin. Phys. C 38, 090001 (2014), 
Sec. 38. 

[4] Numerical Recipes: The Art of Scientific Computing, 
W. H. Press et al, 3rd Edition (2007), Cambridge Uni¬ 
versity Press, Sec. 15.6. 

[51 P. Jimenez-Delgado and E. Reya, Phys. Rev. D 89, 
074049 (2014). 

[6] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, 
Eur. Phys. J. C 63, 189 (2009). 

[7] J. F. Owens, A. Accardi and W. Melnitchouk, Phys. Rev. 
D 87, 094012 (2013). 

[8] J. Gao et al, Phys. Rev. D 89, 033009 (2014). 

[9] V. Radescu, arXiv: 1308.0374. 

[10] S. Alekhin, J. Bliimlein and S. Moch, Phys. Rev. D 86, 
054009 (2012), 

[11] L. W. Whitlow et al, Phys. Lett. B 282, 475 (1992); 
L. W. Whitlow, Ph.D. thesis, Stanford University, SLAC- 
357 (1990). 

[12] F. M. Steffens, M. D. Brown, W. Melnitchouk and 
S. Sanches, Phys. Rev. C 86, 065208 (2012). 

[13] F. Lyonnet et al, JHEP 1507, 141 (2015). 

[14] https : //users . hepf orge . org/~pj imenezdelgado/ 
chromopolis.html. 

[15] M. Arneodo et al, Nucl. Phys. B483, 3 (1997). 

[16] F. D. Aaron et al, JHEP 1001, 109 (2010). 

[17] J. J. Aubert et al, Nucl. Phys. B213, 31 (1983). 

[18] S. J. Brodsky, arXiv:1401.5886 [hep-ph], 

[19] H. Abramowicz et al, Eur. Phys. J. C 73, 2311 (2013). 

[20] S. J. Brodsky, P. Hoyer, C. Peterson and N. Sakai, Phys. 
Lett. B 93, 451 (1980). 

[211 J. Pumplin, H. L. Lai and W.-K. Tung, Phys. Rev. D 75, 
054029 (2007). 

[22] S. Dulat et al, Phys. Rev. D 89, 073004 (2014). 

[23] R. D. Ball et al, Nucl. Phys. B874, 36 (2013). 

[24] T. J. Hobbs, J. T. Londergan and W. Melnitchouk, Phys. 
Rev. D 89, 074008 (2014). 



